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In the scalar tensor gravitational theories Newton's constant Gjv evolves in the expanding universe. 
Likewise, it has been speculated that the acceleration scale oo in Milgrom's modified Newtonian 
dynamics (MOND) is tied to the scale of the cosmos, and must thus evolve. With the advent of 
relativistic implementations of the modified dynamics, one can address the issue of variability of 
the two gravitational "constants" with some confidence. Using TeVeS, the Tensor- Vector- Scalar 
gravitational theory, as an implementation of MOND, we calculate the dependence of Gn and ao 
on the TeVeS parameters and the coeval cosmological value of its scalar field, <^c. We find that 
Gn, when expressed in atomic units, is strictly nonevolving, a result fully consistent with recent 
empirical limits on the variation of Gn- By contrast, we find that ao depends on (pc and may thus 
vary with cosmological epoch. However, for the brand of TeVeS which seems most promising, ao 
variation occurs on a timescale much longer than Hubble's, and should be imperceptible back to 
redshift unity or even beyond it. This is consistent with emergent data on the rotation curves of 
disk galaxies at significants redshifts. 

PACS numbers: 98.80.Es,04.50.Kd,95.35.-|-d,04.80.Cc 
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I. INTRODUCTION 

The debate over the constancy of physical constants 
has simmered ever since Dirac ennunciated the Large 
Numbers hypothesis: very large (or small) dimension- 
less universal constants cannot occur in the basic laws of 
physics P]. In particular, since the dimensionless gravita- 
tional constant is very small, the possibility of variation 
of Newton's constant Gn was raised. The Brans-Dickc 
theory of gravitation among others, can describe such 
variation by adding to the Einstein-Hilbert action an ac- 
tion for a scalar field. 

The past few decades have witnessed extensive searches 
for evidence of variation of some of the fundamental 
constants. Among the methods used have been astro- 
physical observations of the spectra of distant quasars, 
searches for variations of planetary radii and moments 
of inertia, investigations of orbital evolution, searches for 
anomalous luminosities of faint stars, studies of abun- 
dance ratios of radioactive nuclides, and (for current vari- 
ations) laboratory intercomparison of precise clocks 
To cite one example, current data on elemental abun- 
dances, when compared with the theory of Big-Bang nu- 
cleosynthesis, limits the mean rate of variation of Gn 
since early epochs to (Gn/Gn) < 3 x 10~^^yr~^ [3|. 
Obviously any proposed new theory of gravitation must 
be in harmony with this constraint. 

Gat is not the only gravity linked "constant" which 
might be variable. Milgrom's modified Newtonian 
dynamics-MOND f^, which was proposed to explain 
mass discrepancies in galactic dynamics without calling 
on dark matter, introduces a new fundamental parame- 
ter, ao, with dimensions of acceleration. In MOND, New- 
ton's second law a = — V<f>Ar is replaced by 



where $jv is the usual Newtonian potential due to the 
baryonic matter alone, and the function fl{x) smoothly 
interpolates between fl{x) = a; at a; <C 1 and the Newto- 
nian expectation p.{x) = 1 at a; ^ 1. This phenomenolog- 
ical relation, with ao — 10~^°m/s^, has had great success 
in explaining the rotation curves of disk galaxies using 
only the distribution of visible matter, as well as the slope 
and observed tightness of the TuUy-Fisher relation, which 
correlates the luminosity (or baryonic mass) of a disk 
galaxy with its asymptotic rotational velociW J6|l. Recent 
reviews of MOND may be found in Refs. 0, M,M, 

Use of MOND immediately raises a question: are the 
parameters Gn and ao appearing in it constants of na- 
ture, or are they subject to spacetime changes, as is the 
Gn in Brans-Dicke theory? 

Milgrom [5] noticed that the observed value of ao is 
quite close to cHq where Hq is the present epoch Hub- 
ble "constant" . He thus conjectured that Oo may de- 
crease together with the Hubble parameter on cosmo- 
logical timescales @. By contrast Sanders [ll| found 
that in the framework of the tensor-vector-scalar the- 
ory (BSTV) he proposed to provide a cosmological basis 
for MOND, ao grows with time, and would differ signif- 
icantly at redshift z ~ 1 from its present value. This 
would imply significantly reduced asymptotic rotational 
velocity for distant galaxies. Sanders also remarked that 
in BSTV, just as in scalar-tensor theories, Gn varies at 
a rate marginally in conflict with current observational 
bounds on Gn/Gn- 

TeVeS, a new relativistic theory of gravity, was pro- 
posed by one of us [l^] as a basis for MOND. It has been 



explored and subjected to a wide battery of tests 13 1, 
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and has also been extended in various directions 
TeVcS has MOND as its weak potential, low acceleration 
limit, while its weak potential, high acceleration limit 
(1) is the usual Newtonian gravity. TeVeS is endowed with 
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three dynamical gravitational fields: a scalar field (j), a 
timelike unit normalized vector field Ua , and the Einstein 
metric gap on which the additional fields in the theory 
propagate. The theory also employs a "physical" metric 
ga/3 on which gauge, spinor and Higgs fields propagate. 
It is related to gaf3 by 

= e^'^'^gafj - 2uaUp sinh(2(/)). (2) 

The index of Ua or of is always raised with the metric 
5"^, the inverse of gap- 

The equations of motion for the fields in TeVeS derive 
from a five-term action depending on four parameters: 
the fundamental gravity constant G, two dimensionless 
parameters k and K and a fixed length scale ^. We use 
here the form of the action given in Ref. [l^. Variation of 
the action with respect to g""^ yields the TeVeS Einstein 
equations for gap'- 

Gap SttG {fap + (1 - er^'^) u'^ff^(aUp) + T^/j) + Oap- 

(3) 

The sources here are the usual matter energy-momentum 
tensor Tap (related to the variational derivative of Sm 
with respect to g"^), as well as the energy-momentum 
tensors for the scalar and vector fields, 

- ^(2/) / , , ^, , \ ^{y)9aP f.. 

Sap = K (^g^"'U[f,^a]U[u,p]- ^9'^''9'"'U[a,f,]U[r,u]9aP^ , 

- XUaUp (5) 

with 

= ^'(2/); y = kfh'''cb^^cl>,s- (6) 

Each choice of the function ^{1/) defines a separate TeVeS 
theory. Its derivative functions somewhat like the 
fl function in MOND. For y > 0, ~ 1 corresponds 
to the high acceleration, i.e., Newtonian, limit, while the 
limit < fi{y) <^ 1 corresponds to the deep MOND 
regime. We shall only consider functions such that > 
and fi > for either positive or negative arguments. 

The equations of motion for the vector and scalar fields 
are obtained by varying the action with respect to (f) and 
Ua, respectively. We have 

a].p^kG [g''^ + (1 + u'^u''] tap , 

(7) 

for the scalar and 

= SttG (l - e'^^) g'^'^u^f^p - (8) 

for the vector. Additionally, there is the normalization 
condition on the vector field: 

u''Ua=gapu''u^ (9) 



The A in Eq. ([8]), the lagrange multiplier charged with 
the enforcement of the normalization condition, can be 
calculated from the vector equation. 

The three parameters, fc, K and ^, all specific to TeVeS, 
are constant in the framework of the theory, as is G, 
the fundamental gravitational coupling constant, (which 
need not coincide with Newton's Gat). As shown in 
Ref. [H, the measurable quantities Gat and ao can be 
expressed in terms of k, K, i and G. However, that calcu- 
lation of Gat neglected the nonzero cosmological value (pc 
of the scalar field in the scalar equation's matter source, 
thus obtaining spurious cosmological evolution of Gtv- In 
this paper we carry out the calculations in great detail, 
and show that TeVeS predicts a strictly nonvarying Gn 
and only a weak cosmological evolution of ao. This is in 
full agreement with available observational constraints on 
the evolution of Gat, as well as the emerging constraints 
on evolution of ao. As extended rotation curves of high-z 
galaxies become available in the future, they will make 
possible a serious check of the TeVeS's prediction that ao 
evolves weakly. 

In Sec. In] we clarify the sense in which Gn turns out 
to be constant by contrasting physical (atomic) units of 
length with the Einstein units in which the gravitational 
action looks simple. In Sec IIIII we work in the Newto- 
nian (strong acceleration-weak potential) limit of TeVes 
to calculate Gat in terms of the fundamental constant G 
and the TeVeS parameters K and k- Passing to the weak 
acceleration limit of TeVes we calculate in Sec. IIVI ao in 
terms of the TeVeS parameters £, K and k and the cosmo- 
logical value (t>c- In Sec. |V]we estimate the cosmological 
evolution of ao, first naively by assuming that it is tied 
to that of the Hubble parameter which is taken to evolve 
a la GR, and then by setting a bound on the rate of 
(f) evolution from TeVeS's equations. The latter method 
clearly shows that ao evolves slowly on Hubble's scale. 
Sec. IVII summarizes our conclusions. In the Appendix we 
check our methodology by recovering the accepted Gat 
for the Brans-Dicke gravitational theory. 

Greek indeces run over 0, 1, 2, 3 with = t represent- 
ing time; a partial derivative with respect to t is denoted 
by an overdot. We set c to unity everywhere. 



II. DIMENSIONS AND UNITS 

What does it mean to say that Gat is not evolving? Af- 
ter all, a dimensionfull quantity can be caused to be con- 
stant in spacetime by simply choosing the unit in which 
it is measured to have suitable spacetime variation 16]. 
So the only operationally meaningful statement of con- 
stancy of Gat is that some dimensionless combination of 
physical parameters involving Gn does not evolve. Let 
us thus specify such a constant combination in TeVeS. 

Following Dicke's masterful critique [l^, we choose in 
the present section to regard the metric coefficients as 
carrying dimensions of squared length and to think of all 
the coordinates themselves as dimensionless. As men- 
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tioned, in TeVeS the equations for the material fields 
(spinor, gauge and Higgs fields) take their usual form 
when written on the physical metric gai3- In particular, 
we assume that in the stated formulation the dimension- 
less gauge coupling constants and all elementary particle 
masses are constant in spacetime. What is being assumed 
is that the system of units reflected by gap uses a par- 
ticle mass, say the proton's rup, as its local mass unit. 
Likewise, the physical parameters c and /i, which appear 
in the various equations, are supposed constant. Since 
?i, c and rUp all bear different dimensions, this last re- 
quirement has fixed the system of units (up to the trivial 
freedom to double the length unit everywhere, etc.). The 
statement that Gm is nonevolving in such "atomic" units 
is thus equivalent to the statement that GNirip/hc is a 
spacetime constant. 

The above units differ from those carried by the Ein- 
stein metric gap which is the one used in formulating 
the TeVeS equations. (Here it may prove conceptually 
useful to regard the Tap, which is calculated by varying 
the matter action with respect to g'^^ , as reexpressed in 
terms of gap)- According to Eq. ([2]), for like coordinate 
increments, the physical distance in the space orthogo- 
nal to u" (the space whose metric is gap + UaUp) is a 
factor e"*^ times the distance paced out by the Einstein 
metric itself. By contrast, the physical distance coUinear 
with the timelike vector field is e"^ times that given 
by the Einstein metric. In other words, the Einstein unit 
of length is not only spacetime varying but also space- 
time anisotropic with respect to that in "atomic" units. 
In Einstein units we may still regard c and h (but not mp 
or its corresponding Compton length) as constants. Ac- 
cordingly, we have set c = 1 everywhere. As mentioned 
G is constant in Einstein units. Were we to set it as well 
as h to unity, we would be in Planck units. However, we 
shall refrain from this last step and continue to exhibit G 
explicitly. The main question we shall be asking is, how 
do Gn or ao, when appropriately calculated in physical 
units, relate to the TeVeS constants G, etc.? 

III. NEWTON'S Gn IS CONSTANT 

We begin by showing the relation between Newton's 
constant and the TeVeS coupling constant G. By defini- 
tion, Newton's constant enters through the relation 

$^ = -G^ [ -f^-d^i' (10) 
J |x' - x| 

with m = J p{Sc')d^x' the physical mass and x the Carte- 
sian coordinate that marks physical distance. Since we 
expect TeVeS to have the customary weak field limit, the 
Newtonian potential (jlOp should enter in the customary 
way in the linearized form of the physical metric (the 
metric measured by instruments made of matter). We 
thus expect that 

d~S^ = -(1 + 2$Ar)dT2 + (1 - 2$Ar)di • di. (11) 



where r is the coordinate that marks physical time. 

To be able to compare Einstein and physical metrics, 
we must use the same coordinates for both. To maintain 
consistency with previous work p^ . [isj we choose new 
coordinates t and x in terms of which the asymptotic 
physical metric, though flat, differs slightly from stan- 
dard Minkowsky form. The relations between physical 
distance and time, x and r, and the coordinates x and t 
are 

x = e-"^'=x; T = e'^'=t. (12) 

With this in mind, we can rewrite the weak field limit of 
the physical metric outside its source in the form 

ds^ = -(l + 2$jv)e2'^<=dt2^(l-2$jv)e"2"^'^dx-dx. (13) 

To relate Gjv to G we must solve TeVeS equations and 
use the solutions to construct the physical metric, which 
should turn out to be identical to Eq. Consider the 

following ansatz for the Einstein metric to first order: 

ds'^ ^ -{l + 2V)dt^ + {l~2V)dyi- dy.. (14) 

The vector field must be timelike; we shall take our 
coordinate system to coincide with the "rest frame" es- 
tablished by u". In view of the normalization condition 
(O, the vector field is given to first order by 

"a = {-(l + ^^),0,0,0}. (15) 

The scalar field may be written 

= 0, + 5dp, (16) 

with (pc the nonzero cosmological value of cf) and (50 rep- 
resenting the local departure from it; 50 is to be regarded 
as of the same order of smallness as V . 

Now for the details. To find the potential V , we need 
to solve the Gtt Einstein equation to first order. We take 
for the energy-momentum tensor the familiar ideal fluid 
form 

Tap = pVaVp + p {gap + VaVp), (17) 

where p is the proper energy density, p the pressure and 
Va the 4-velocity, all three referred to the physical metric. 
We assume that the fluid is stationary in the coordinates 
chosen, so the spatial part of v" must vanish, i.e., Va 
must be parallel to Ua- When account is taken of Va 
normalization with reference to gap we get (T2| 

Va = e'^Ua. (18) 

In the nonrelativistic approximation p is negligible rela- 
tive to p; thus 

Tap = pe^^UaUp. (19) 

We now substitute ([llD , HH]) , (HH), and into the 
TeVeS equations, and retain only first order in V and 
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5(j)- We start by solving the temporal component of the 
vector equation ([5]) for A (the other components are zero 
to linear order), obtaining 



A = -KW^V - 167rGpsinh(2(/)c). 



(20) 



When this is substituted in the Gu Einstein equation, it 
becomes, to first order, 



AtiG 
l-K/2 



(21) 



We have neglected in the r.h.s. terms of order 5(j) which 
would source terms of second order in V . With the 
boundary condition y — > for |x| — > oo the solution 
is 



V = 



a 



1 - KI2 



(22) 



in which integral p is regarded as a function of the coor- 
dinates X, not of the physical distances. 

This result must make its way into the physical metric. 
Using transformation ([2]) in Eq. ([14]) we get 

ds'^ = -(1 + 2V)e^'>'df + (1 - 2y)e"2'*dx • dx, (23) 

which becomes-to lowest order in (50- 

d§^ ^ -{l + 2V+2S(j))e'^'l'^dt^ + {l-2V-25(t))e-'^'''-'dx-dx. 

(24) 

Comparing with Eq. (|13p we may identify the Newtonian 
potential 



N 



V 



(25) 



We now need only calculate d(j>. 

In the scalar equation ([7]) we substitute ^ = 1 because 
we are concerned with the Newtonian limit in which fi 
approaches unity; any small corrections to it may be dis- 
carded since we work here to first order in Scb: 



(26) 



Note the factor e ^'^'=, the correct asymptotic value of 
e~ 'f' . which was missed in Eq.(53) of Ref. The solu- 
tion of this last equation in accordance with the boundary 
condition (h ~^ (hr as x. oo is 



^^d^x'. 



(27) 



47r J |x' — x 
Substituting Eqs. (HH) and (QZl) into Eq. (gS]) we get 



G 



^d'x'. 



(28) 

Finally, we use relation to switch back to physical 
length coordinates x: 



$Ar(i) 



where p(x) = p{xe~'^'') is the energy density distribution 
in physical units. Comparing with Eq. (jlOp we obtain 
Gjv in terms of the TeVeS parameters: 



' {2-K)k ' 



(30) 



Thus the ratio G/Gn turns out not to depend on 0c, 
the asymptotic cosmological value of 0, a cosmologically 
evolving quantity. Since G, k and K are constant param- 
eters, Gn (as measured in physical units) does not evolve 
with cosmological epoch. Further, since Gat is observa- 
tionally positive, while it is natural to expect that G > 0, 
we must restrict K to the ranges K<2orK>2 + 87T/k. 
It is amusing that one can reconcile the small observed 
value of Gat (more properly of GNTn^/K) with strong 
gravity {Gm^/h which is not especially small) for the 
family of TeVeS theories for which K is only very slightly 
above the critical value 2 -|- 87r/fc. In such a theory the 
true Planck length (G/i)^/^ could be commensurate with 
elementary particle scales like h/mp or be even larger, all 
this without resorting to brane physics. 

The result that Gat is constant is surprising in view of 
the fact that TeVeS contains a scalar sector. To check 
our methodology we work out, in Appendix A, Gat for 
Brans-Dicke theory by following the track set out in the 
present section. We obtain the accepted law of evolution. 



IV. THE MOND ACCELERATION SCALE 

The MOND acceleration scale ao can also be calculated 
in terms of the TeVeS parameters. Since MOND is the 
small acceleration, weak potential limit of TeVeS p^ . 
we can again work with the linear approximation to the 
physical and Einstein metrics; however, in the present 
case /X < 1. As in Ref. [H we shall look for a MOND- 
like equation of the form ^ . We shall then attempt to 
identify the MOND function jl and the combination of 
TeVeS coupling constants which is equivalent to ao. For 
simplicity, we shall assume spherical symmetry; it can be 
shown that our results hold for asymmetric systems as 
well [12. 

We write the Einstein metric for weak potentials ex- 
actly as in Eq. ([T^ . while in contrast to Eq. ([TH)) the 
physical metric is expected to be 

d~s^ = -(1 + 2$)e2'^'=di2 + (1 - 2$)e-2<^'=dx • dx, (31) 

with <I>, the MOND gravitational potential, replacing $ at . 
By transforming from Einstein to physical metric in ac- 
cordance with Eqs. ^ and (fT5 |) - pB|) we find to first order 
that 



^ = V + 5(t). 



(32) 



The contribution of y to $ is the same as that to 
because the terms depending on jj, in the Gu equation 
from which V arises are all of second order in 56, and 
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thus stand for higher order corrections. Thus we have V 
as in Eq. 

In determining Scj) here we must take into account the 
fact that in the weak acceleration hmit /i < 1. The scalar 
equation ([7]) takes the form 



V • 



(33) 



Comparing this equation with Poisson's and using Gauss' 
theorem in the spherically symmetric case gives 



V(50 



47r/i 



(34) 



Then in view of Eqs. and ^) the gradient of the 
total potential $ satisfies 



1 k \ 

\~K/2 ^4^) 



(35) 
(36) 



Eq.dnnD is the desired MOND-hke equation. To find 
Milgrom's parameter ao we proceed to the extreme 
MOND regime defined by the condition /i <C fc/(47r). 
There Eq. p6p gives jl ~ A'K^/k. Substituting this in 
Eq. ([55l) and comparing the result with Eq. reveals 
that in the said limit V$ ~ V^0. This implies that in 
Eq. (gSl) W is then negligible. 

Instead of focusing on the toy form of /i(y) from 
Ref. [H, or any other ansatz for it, let us be very gen- 
eral. We just require that for < y <C 1, ^Ji{y) « D^Jy 
where D is some positive constant. Going to sufficiently 
small y so that /j, ^ A:/(47r) we have, in view of the last 
paragraph, that 

jl « 47r^/fc « A-KDk^^''^l\V64>\ « 47ri:)fc~^/2£|V$|. 

(37) 

Substituting this in Eq. ([HS)) and transforming all ocur- 
rences of x (including in the gradients) to x by means of 
Eq. ini), we get the extreme MOND equation 



with 



Oo = 



G ^ke'l"^ 
Gn 47rD^ 



(38) 



(39) 



We have here employed the definition (dU]); it is under- 
stood that <I> also is to be regarded as a function of the 
physical length coordinates x. We see that Milgrom's 
acceleration scale ao depends on the TeVeS parameters 
k, K and I (all constant in Einstein units) as well as on 
the constant coefficient D associated with the function 
J-{y). But unhke Gjv, is predicted in TeVeS to evolve 
cosmologically in consonance with e'^'= . How fast an evo- 
lution it is capable of is the subject of the next section. 



V. PREDICTED EVOLUTION OF oo 
A. The naive MOND viewpoint 

The pure MOND paradigm is arnbiguous about the 
time evolution of Oq. Milgrom llCll | remarks on the 
numerical coincidence between ao and the observed c_ffo 
{Hq is the present value of the Hubble parameter H) or 
Oo and the value of the cosmological constant A inferred 
from the acceleration of the cosmos. If the former coinci- 
dence bespeaks of a physical connection, then one would 
expect cosmological evolution of ao with ao oc iJ, while 
if it is the second coincidence that properly reflects the 
physics, then ao should be strictly constant. 

How big an evolution to expect in the first case ? Since 
naive MOND does not provide a consistent cosmology, 
we shall here use cosmology a la GR. Let us write the 
Friedmann equation in GR for a cosmological model with 
curvature index k: 
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A 
3" 



(40) 



Here 6 = b{t) is the expansion factor with value 6o at the 
present time, pmo is the present value of the mass density 
of pressureless matter, and we are neglecting radiation's 
contribution because we focus on the more recent uni- 
verse. Differentiating with respect to t and dividing out 
by 2H gives 



H__ _( ATlGprnghl 



H. 



(41) 



As customary we may introduce densities as fractions 
of the present critical density 



SnGpmo 



A 

3^' 



(42) 



so that Q.ra~\-^K.+^K = 1 ou account of Eq. (I^O]) evaluated 
at the present epoch (when 6 = 60 and H — Hq). We see 
that 



iH/H)c, 



(43) 



The standard cosmological model obtains values for the 
ri's from various observations, e.g., those of the cosmo- 
logical microwave background anisotropy spectrum. 
comes out either zero (flat space) or positive (hyperbolic 
space) and very small on scale unity. By contrast ri,„, 
which includes the contribution from putative dark mat- 
ter, is assigned a value of about 0.25. We may thus 
conclude that at present do/ao, which is the same as 
{H/H)o, should be about — 0.25i?o- Thus the present 
day timescale of ao variation is four times longer than 
the Hubble scale. 

As we go back in time ao should scale proportionately 
to the coeval H. We may recast Eq. (|40l) as 



H ^Ho + z)2 + + n„,ii + zf 



1/2 



(44) 
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with 1 + z = bo/b. With « 0, « 0.25 and 
~ 0.75 as in the standard model, the curvature term 
in the square brackets in the last equation remains neg- 
ligible, while by z « 1 the matter term will have come to 
dominate the A term. We then have for z > 1 



ao(z) « a(0)(l + z) 



3/2 



(45) 



which implies a drastic change of ao between z of a few 
and today. 

However, it could be claimed that to keep in the spirit 
of the MOND paradigm one should, apart from retaining 
Hk, ~ 0, equate ^Im with the baryon fraction fib = 0.04 
inferred in standard cosmology. This last can easily ac- 
commodate still unobserved massive neutrino or baryonic 
matter which is nowadays invoked in MOND in connec- 
tion with the large clusters of galaxies d, [13, [3. Of 
course to be consistent we should then put JIa ~ 0.95. 
With this set up the matter term in Eq. becomes 
comparable with the A term only for 2^2, and ao will 
follow the law (|l5l) for z > 2. For z <C 1 we would have 
from Eq. ([^5]) that Oo changes on a timescale 16 times 
longer than Hq^. 

The above discussion is instructive; but it is hardly 
trustworthy as underlined by the contrasting results it 
can yield. The crux of the problem is, of course, that 
MOND is not a nonrelativistic limit of GR, yet this last is 
being used to work out the cosmology. This inconsistency 
can be avoided by calculating the cosmological evolution 
of Oo entirely within TeVeS, which does have MOND as 
a nonrelativistic limit. 



B. The TeVeS viewpoint 

We found in Eq. ([55)1 that ao, as defined by small scale 
MOND dynamics, has a e'^ dependence, where by 4> is 
meant 0's cosmological value 4>c- We are thus invited 
to establish the cosmological evolution of cj). It will be 
useful to distinguish here, as we did in Sec. Illli between 
the coordinate time t and the physical time t. 

The Einstein metric for a Friedmann-Robertson- 
Walker model is 



sm 



}dip^)] (46) 



where /(x) is either x (open model with flat spaces, 
K = 0) or sinhx (open model with hyperbolic spaces, 
K — —1). As in Sec. VII of Ref. [13 we shall take 
= {1,0,0,0} and (j) = (pit), consistent with the time- 
like character of the vector and the assumed isotropy and 
homogeneity of space. Then according to Eq. ^ we ob- 
tain the physical line element ds^ by multiplying the tem- 
poral part of gap by e^"^ and the spatial parts by 



,-20. 



d~s^ = -dT^ + b{tf[dx^ + fixfide^ + sin^ed^^)], (47) 
with 

b = e-'^b- dr = e'^dt. (48) 



The r here is the physical time since it acts as proper 
time of commoving observers, c.f. Eq. (H^. 
From ao cx; e*^ it follows that 



c?ao /dr d(j) 



The first integral of Eq. ([7]) for 
the case of ideal fluid matter: 



(49) 



is given in Ref. [13 for 



^(-2fcr0^)(/) 



-k 

2P 



G{p + ■ip)e-^H^dt, (50) 



Here, as earlier, an overdot designates a derivative with 
respect to i, not r. Since the physical energy density, 
p, the physical pressure, p, and the TeVeS parameter 
k are all positive, we see that cj), and consequently also 
d(j)/dT, are negative. Thus by Eqs. ao is strictly 

decreasing with physical time t. But because the integral 
above includes contributions from early times when p is 
not negligible, and because of the complicated factor /x, 
the said equation is far from convenient for estimating (p. 
We shall instead estimate d(t)/dT by way of the Einstein 
equations. 

First we compute the physical Hubble parameter H: 



H 



db/dr 



dr' 



(51) 



Next we compute A from the vector equation ([5]); it takes 
the form [11 



A = Stt [^i^Vfc - 2G sinh(20)/5] . 



(52) 



Finally we write down Einstein's equations ^ sans the 
cosmological constant and with a perfect fluid as matter. 



62 

62 



K, 
62 



-pe 



-24> 



3fc 3fc2£2 ' 



(53) 



where y here is identical to the argument of p in Eq. (j50p . 
Taking into account that J- > and p > and that 
(f) < 0, we see that for a spatially flat or hyperbolic cos- 
mological model (k < 0) 



> 



/I67r/x\ 



1/2 



dr' 



(54) 



This could be a strong inequality if the p term in Eq. 
is dominated by the matter energy density. In any case 
from Eq. (jST]) we see that 



dr 



< H. 



(55) 



It is clear from Eqs. ((49|) and (|55p that for any choice 
of p-so long as it is positive. 



dao / dr 



< H 



db/dr 



(56) 
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Thus within any reasonable TeVeS theory ao's evolution 
is slower than the Hubble expansion at the same epoch, 
and it can be much slower, provided only fi is not small 
compared to unity in recent epochs and A: < 1. 

A case in point is the TeVeS theory investigated in 
detail in Ref. 12. It incorporates a function J-{y) for 
which > 1 for y < 0. As shown there one then 

needs fc ^ 1 for TeVeS cosmology to be consistent with 
causality. Thus by Eq. dSS]) \d(j)/dT\ < H. Then by 
virtue of Eq. this implies 



dao /dr 



do 



db/dT 



(57) 



Thus at all epochs the evolution of Qq occurs on a 
timescale much longer than Hubble's. Put another way, 
as one goes back in time, ao{z) grows much slower than 
bo/b{z) or 1 + z. 



VI. CONCLUSIONS 

In this work we have calculated Newton's constant 
and the MOND acceleration scale ao in terms of TeVeS' 
parameters. We find that Gn does not depend on the 
dynamical scalar field of the theory, and is thus strictly 
constant in cosmology. This corrects an impression that 
one might obtain from Ref. It also shows that analo- 
gies drawn between TeVeS and familiar scalar-tensor the- 
ories can lead to incorrect inferences. Our result agrees 
with known facts: all existing data point to a nonvarying 
Gn i. 

We also find here that in a cosmological setting ao 
varies as the exponential of the scalar field, thus de- 
creasing with time. However, a detailed consideration 
of TeVeS isotropic cosmological models strongly suggests 
that the Qq variation occurs on scales much longer than 
the Hubble scale. This result is in contrast to a naive view 
which regards Oq as physically connected to the Hubble 
parameter; in such eventuality ao variation would most 
likely occur on the Hubble scale (we have discussed in- 
evitable ambiguities in this point of view). 

At present there are not enough quality data to test 
the TeVeS prediction of slow ao evolution. Clues as to 
the evolution of ao could be gleaned from existing data 
on the TuUy-Fishcr relation at epoch z ~ 1. The TuUy- 
Fisher relation in the form = GnUoM, with M the 
total baryonic mass of the galaxy and Voo its asymptotic 
rotation velocity, emerges naturally in MOND. Evolu- 
tion of ao would entail evolution of the coefficient in the 
TuUy-Fisher relation or, equivalently, of the zero point of 
the plot of logM vs logUoo for disk galaxies, The mea- 
ger available data are consistent with no evolution of the 
TuUy-Fisher relation back to z « 0.6 [l^. In addition, 
Milgrom's MOND analysis [13] of recent data by Gen- 
zel, Tacconi et al. [l^l on the rotation curve of a galaxy 
at z — 2.38 seems to be consistent with unchanging ao 
(although that rotation curve does not extend as far as 
would be desired for this kind of an inference). 
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APPENDIX A: CALCULATION OF Gn IN 
BRANS-DICKE THEORY 

We show here that the methodology of Sec. IIIII will 
yield familiar results when applied to a pure scalar-tensor 
theory such as Brans-Dickc (BD) theory. Whenever fea- 
sible we shall couch the equations in the notation of 
Sec. mi 

Following Dicke [l^ we transform the BD gravitational 
action [3] to the Einstein frame; we shall, however, leave 
the matter action in the physical frame in parallel with 
our treatment of TeVeS [ij] : 



S = 



1 



167rG' 



i^-i(2. + 3)^ 



{-g)'/^d^x 



+ C^r^i-gy^'dS 



(Al) 



In the above co is the celebrated BD parameter and A, 
a dimensionless entity, represents the BD field in units 
of the fundamental constant G~^, i.e., A = Gcj). The 
first line of the action is stated in terms of g^fs while the 
matter action takes its usual form when written in the 
gap metric. In BD theory 



Jaf3 — A ^gap; 



(A2) 

On account of the definition of the matter's energy mo- 
mentum tensor as a variational derivative of the matter 
action we have 

-2 5[C^{^~gf/^] = f^p{~~gf/^5r^ 
= {-gf'^[fap5g'^f'\~^ + fapr^\-H\l (A3) 

where the second line results on account of the transfor- 
mations ()A2|) . Now variation of gap in S together with 
the identity 

5[R{-gf'^] = Gap{~gY'^5g''l^ +honnda.ry terms (A4) 
and our last result yields the gravitational equations 

Gap — (2t^ + 3)A^^ ^,a^,P ^,p.^'^9ap 



SirGTapX ^, 



(A5) 



the counterpart of our Eqs. ([3]), while variation with re- 
spect to In A yields the BD scalar equation in the form 



^[,"^(lnA),„(-,)V^],,= 2^ 



(A6) 
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which is the counterpart of Eq. ([7]) . 

Let us solve the equations for a stationary situation to 
hnear order by writing 



9af3 = VaP + hai3 and A = Ac + C, 



(A7) 



with r]ap and Ac the asymptotic values of the Einstein 
metric and scalar field (where spacetime is assumed 
flat). Now according to Eq. (fT7|) for ideal fluid matter 
Taf3 g"^ = —p + Sfj. In first approximation we may ne- 
glect the p. Then to first order in ^ and (and neglect- 
ing any temporal variation of the cosmological boundary 
value Ac) Eq. (|A6p takes the form 



SttG/Ac 



2uj + 3 

whence in analogy with Eq. 



c 



2G/Ac f /5(x') 



2u; 



X . 



(A8) 



(A9) 



This last result shows that the first term in the r.h.s. 
of Eqs. (jASp is of second order in Gp, and thus negligible 
compared to the matter term. Again, from Eq. (|17p we 
see that here Taf3 ~ pVaVp. If the matter is static we 
have from the normalization of in the physical frame 
that = —gtt^a^p- This will also be true to good 

approximation if the matter flows in space provided that 
Wet's spatial part v is small compared to unity (errors will 
be of ©(v^)). Hence the BD gravitational equations are 

G^^ « -SnGX-^pgu S'Jj^ « SirGX-^pS^Sj^, (AlO) 



where we have used relations (jA2l) and dropped from the 
last expression, subdominant terms with extra factors 
of ha/3 and ^. These equations are just the GR Einstein 
equations in metric gaf3 for a quasistatic mass-energy dis- 
tribution p, but with G/X~^ playing the role of gravita- 
tional constant. 

We know that to flrst order such Einstein equations 
have the line element (I14p as solution with V signifying 
the usual Newtonian potential. In light of our remark 
about the gravity constant we must write here the fol- 
lowing analog of Eq. (|22|) : 



V = -GA; 



P(x') 



(All) 



Using the transformations (|A2|) we evidently have to first 
order in V and that 

ds^ = -il + 2V - C/Xc)Xc^dt'^ 

+ (l-2F-C/Ac)A^Mx-dx (A12) 



In order that the physical line element be asymptoti- 
cally Minkowski, we must define, in analogy with rela- 
tions (I12p . the physical time r and physical lenght coor- 
dinates x: 



X = A;1/2x; r = X-^/H. (A13) 
The line element here has thus a form that contrasts that 
of Eq. HH) for GR and TeVeS: 



ds 



Here 



-(l + 2$w)dr2 + (l-2K$Ar)dx-di. (A14) 



Gn 

K 



V--C/X. 

G 2cj + 4 
aI 2w 3 ' 

UJ + 1 



-Gn / ^^d'i', (A15) 
x' - X 



(A16) 
(A17) 



In Eq. (jAlsp we have absorbed one factor Ac ^ into the 
integral to convert from x to x. 

Comparison with Eqs. piI)l - (fTT|) shows that Gat here 
is properly regarded as the Newtonian gravity constant. 
Our Gn concurs with Brans and Dicke's [4I showing 
clearly that in BD theory the Newtonian gravity "con- 
stant" , by virtue of its strong A dependence, evolves cos- 
mologically, in contrast to the case of GR or of TeVeS. 
Our value for the coefficient k also coincides with that 
obtained by Brans and Dicke [3]. The fact that k ^ 1 is 
responsible for gravitational lensing being smaller in BD 
theory than in GR, is again in contrast to the TeVeS case 
for which gravitational lensing is the same as that in GR 
for the same source p(x). 
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